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The relevance of zero-energy functions, coming from zero-energy modes and present 
in the structure of bosonic Green's functions, is often underestimated. Usually, their 
values are fixed by assuming the ergodicity of the dynamics, but it can be shown that 
this is not always correct. As the zero-energy functions are connected to fundamental 
response properties of the system under analysis (specific heat, compressibility, suscepti- 
bility, etc.), their correct determination is not an irrelevant issue. In this paper we present 
some results regarding the zero-energy functions for the Heisenberg chain of spin- 1/2 with 
periodic boundary conditions as functions of the number of sites, temperature and mag- 
netic field. Calculations are pursued for finite chains, using equations of motion, exact 
diagonalization and Lanczos technique, and the extrapolation to thermodynamic limit is 
studied. 

Introduction The issue of ergodicity of physical systems is not new. A system is ergodic 
if it goes through every point of phase space during its time evolution. In such systems, the 
equilibrium averages (i.e., the time averages) are equal to the ensemble averages, which are 
much easier to compute. Anyway, we have to face the problem of non-ergodicity of many 
physical systems (e.g., the existence of even one integral of motion divides the phase space into 
separate subspaces not connected by the dynamics) . The lack of ergodicity has measurable 
effects as the well-known difference between the static isolated (or Kubo [1]) susceptibility 
and the isothermal one [2, 3]. 

In the widely used formalism of Green's functions the issue of ergodicity appears as a diffi- 
culty in the determination of the zero-frequency functions present in the bosonic propagators 
[4]-[ll]. The problem reappears also in e.g., composite operator method [12]. The equations 
of motion do not uniquely determine causal Green's functions and correlation functions but 
only up to some momentum function which severely affects the self-consistent calculation of 
the retarded Green's functions too. Usually, these zero-frequency functions are fixed by as- 
signing them their ergodic values, but this can not be justified a priori. Wrong determination 
of them dramatically affects the values of directly measurable quantities like compressibility, 
specific heat, magnetic susceptibility. According to this, they should be calculated case by 
case. 

In this paper, we calculate the zero- frequency functions for the Heisenberg chain of spin- 1/2 
and show that they take non-ergodic values for finite lengths and in the bulk limit too. 

Definitions One of the main aims of condensed matter calculations is to compute the 
correlation functions of the physical system under analysis: 
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where i = (i, t) stands for both spatial i and temporal t coordinates and tp for any operator 

relevant to the dynamics under study. The Fourier transform in frequency of the correlation 
functions can be expressed in terms of the eigenstates |n) and eigenenergies En of the system 
in the following way 

C{i,j,uj) = 27r,5Hr^^t(i,j) 

+ ^ <i~'"'"{nmi)\m}{m\i^^mn)SliJ + {En-Era)] (2) 

where the zero-frequency function r^y,t(i,j) is defined as: 

r^^t(ij) = | Yl e-''^"HV(i)|m>(m|V't(j)|n) (3) 

n,m 

Its ergodic value is [1, 13]: 

r;7t(i,j) = (^(i)>(^t(j)) (4) 

In this manuscript, we focus on the Heisenberg chain of spin-1/2 with periodic boundary 
conditions in presence of an external magnetic field directed along z and proportional to h: 

H = jJ2sii)S{i")-hJ2s^i) (5) 

i i 

whore i" = (i+l,t), J is the exchange coupling and S{i) is the spin-1/2 operator at the site 
i of the chain of N sites. 



Results We diagonalizc exactly the Heisenberg model on three sites in magnetic field, ob- 
taining eigenvalues and eigenvectors. Then, Eq. (3) yields the following functional dependence 
on temperature and magnetic field of the zero-frequency function rs^s^(i,j) (see Fig. 1): 



Ts^s, (i,i) 



rs,s,(i,i + l) = 



I_2e'5''-Ke2'"' 



36 



36 



+ 



5 ^ 4e/3'' + Se^'^'' 



18 



l-he2/3''-K2ei'3''e'3ft 



(6) 



(7) 



We can see that the zero-frequency functions behave quite differently then their ergodic values 
M'^ = {SzY £^nd with a quite strong dependence on both magnetic field and temperature. 
According to this, a careful determination of them is absolutely necessary and we cannot rely 
at all on the ergodic values. 

Let us now analyze the situation as function of the number N of sites. To perform calcu- 
lations for chains with a number of sites greater than 4 it is necessary to use numerical tools 
like the exact diagonalization and the Lanczos method. The thermodynamic limit will be 
extrapolated by means of a 1/N analysis. 

For zero temperature the ergodic value of the on-site zero- frequency functions Tg+g- is 
always zero, not depending on the magnetic field; for the on-site Vs^'S'^ it is zero excluding 
three cases: 1) J < and finite magnetic field (F^'^^g^ = i); 2) J > and high magnetic field 
h ^ J (rif/^^ = i); 3) J > and magnetic field with values corresponding to energy level 
crossing (for finite N; e.g., h = 3 J/2 for 3 sites). 
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Figure 1: For the 3-sitc chain, the on-site and nearest-neighbor zero frequency functions, 
for antiferromagnetic (,/ > 0; sohd hne and dotted hne, respectively) and ferromagnetic 
(J < 0; dashed and dot-dashed hne, respectively) exchange coupling, together with squared 
magnetization M'^ (double dot-dashed line), at magnetic field h/\J\ = 1, are plotted versus 
temperature 



Properly calculated zero- frequency functions Fs^s^ and Fg+g- behave differently. For 
J < 0, T = and h = both constants are non-zero; this is due to the degeneracy of the 
ground state coming from the fact that both and arc constants of motion and higher 
value of the total spin is selected. We also have the exact relation Fgis^ — ^Tg+g-, which 
can be derived from the algebra in absence of magnetization. The inclusion of a nou-zcro 
magnetic field removes the degeneracy, by fixing the value of and giving dynamics to , 
and switches the functions to their ergodic values. 

For J > 0, T = and h = both constants arc zero or non-zero, depending whether 
the number of sites in the chain is even or odd, respectively. The relation Fazs^ = iF^+g- 
still holds as the magnetization is zero in this case too. The discrepancy with respect to the 
ergodic values for odd chains is due to the impossibility to pair up all spins on selecting the 
lower value of the total spin: even chains have a unique S"^ = singlet-like ground state, the 
odd ones have a degenerate = 1/2 ground state. The inclusion of the magnetic field sets 
Tg+g- again to zero (except for the values generating a level crossing and, consequently, a 
degeneracy), while Tgzgz keeps its previous value or jumps to a higher one passing through 
level crossing values. As the magnetic field becomes high enough the system becomes fully 
polarized, the ground state is non-degenerate as for J < and h > 0, and both constants set 
at their ergodic values. 

These results were obtained solving numerically chains up to 25 sites, using exact diago- 
nalization and Lanczos method. Now the main question: do these results hold also in the 
thermodynamic limit? In Fig. 2 we show the zero- frequency function Tg^s^ as a function of 
the inverse of the number of sites 1/N for J > and J < 0, T = and h — 0. For J < 0, 
Tg^g^ (i, j) is independent of the distance |i— j| and is exactly given by 1/12-1- 1/(6A'^) [14]. For 
J > 0, Fszs2(i, j) is a function of the distance |i— j| and scales with A'^ as shown in the figure 
(the lines represent very accurate fits). According to this, we have a non-ergodic behavior for 
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Figurc 2: Zcro-frcqucncy functions for a Af-sitc chain with ferromagnetic (J < 0; O) 
and antiferromagnetic (J > 0) couphng. In the latter case, zero-frequency functions are 
calculated on-site (A), for first (v) and second neighbors (<0>). The lines are fits of the 
numerical data, according to the formulae displayed in the legend 



J < which remains in the thermodynamic hmit too. We wish to point out that this result is 
exact and must be taken into consideration by whoever needs to compute physical quantities 
like spin correlation functions and susceptibilities in the ID Hcisenberg model by means of 
analytical methods. Fs^s* should be set to 1/12; any other value, the null ergodic one too, 
will induce wrong results and any approximation scheme not taking this into account, directly 
or self-consistently, will surely fail to reproduce the physical properties of the system. Non- 
ergodic behaviors for finite values of temperature and magnetic field cannot be excluded at all 
and are instead probable; we are performing numerical calculations in this direction and the 
results will be published elsewhere. Finally, it is worth noting that non-ergodic behaviors with 
a functional structure (i.e., zero- frequency functions with dependencies on temperature and 
magnetic fields as rich as those found for finite size chains) require the maximum attention in 
building up reliable conserving and/or self-consistent approximation scheme as any mistake 
in reproducing them will mine the possibility to comprehend the low-energy properties of the 
system. 
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